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ABSTRACT 



We present new photometry of 16 local Seyferts including 6 Compton-thick sources in A'-band filters around 12-yum, obtained with 
the VISIR instrument on the 8-m Very Large Telescope. The near-diffraction-limited imaging provides the least-contaminated core 
fluxes for these sources to date. Augmenting these with our previous observations and with published intrinsic X-ray fluxes, we form 
a total sample of 42 sources for which we find a strong mid-infrared:X-ray (12.3 ^m:2-10 keV) luminosity correlation. Performing a 
physically-motivated sub-selection of sources in which the Seyfert torus is likely to be best-resolved results in the correlation Lmir^ 
Lx 1,11±0 i with a reduction of the scatter in luminosities as compared to the full sample. Consideration of systematics suggests 
a range of 1.02-1.21 for the correlation slope. The mean 2-keV:12.3-yum spectral index (»ix) is found to be -1.10±0.01, largely 
independent of luminosity. Observed 12-yum bolometric corrections range over as 10-30 if a known luminosity trend of intrinsic X-ray 
bolometric corrections is assumed. Comparison with ISO data spanning a larger luminosity range suggests that our correlation can be 
extended into the quasar regime. The fact that unobscured, obscured and Compton-thick sources all closely follow the same luminosity 
correlation has important implications for the structures of Seyfert cores. The typical resolution-limit of our imaging corresponds to 
~ 70 pc at a median z = 0.01, and we use the tightness of the correlation to place constraints on the dominance of any residual emission 
sources within these physical scales; we infer an upper-limit of «;40% of the unresolved flux for any contaminating star-formation, on 
average. We suggest that uncontaminated mid-IR continuum imaging of AGN is an accurate proxy for their intrinsic power. 

Key words, galaxies: active - infrared: galaxies - X-rays: galaxies 



1. Introduction 

Typical broad-band spectra of active galactic nuclei (AGN) are 
known to show a significant bump in the infrared regime (e.g. 
Barvainis 1987, Sanders et al. 1989, Elvis et al. 1994, and refer- 
ences therein). This is thought to be due to absorption of intrinsic 
AGN ultraviolet, X-ray and optical radiation by dusty clouds in 
the AGN 'torus' on ~pc-scales, which subsequently re-radiate 
this energy at infrared (IR) frequencies. But isolating the AGN 
emission from other ubiquitous contaminating sources, such as 
bright starbursts within the inner few kpc of AGN host galax- 
ies (e.g. Rodriguez Espinosa et al. 1987; Norman & Scoville 
1988; Gonzalez Delgado et al. 2001), is difficult due to the fi- 
nite resolving powers of IR observatories. Telescopes with im- 
proving sensitivity, wavelength coverage and resolving power 
have changed the situation steadily over the past two decades 
(e.g., Sanders et al. 1989; Mulchaey et al. 1994; Spinoglio et al. 
1995; Krabbe et al. 2001; Lutz et al. 2004; Buchanan et al. 2006; 
Ramos Almeida et al. 2008; Mushotzky et al. 2008; Melendez 
et al. 2008). These works find increasingly-stronger correlations 
between the measured IR powers with the intrinsic (X-ray, bolo- 
metric or various forbidden emission line) AGN luminosities, 
but still with significant scatter. How much of this scatter is in- 
herently due to the physical conditions of the AGN and torus 
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clouds, as compared to observational selection effects, remains 
an important unresolved issue. 

In two papers (Horst et al. 2006, 2008, hereafter Paper I 
and Paper II, respectively), some of us presented high-resolution 
imaging of nearby AGN with the VLT/VISIR mid-infrared cam- 
era (Lagage et al. 2004) in filters covering the 8-13 pm atmo- 
spheric transmission window. The diffraction limit of the 8.2m- 
diameter mirrors of the VLT over this V-band window is a;0''3- 
0"4, which is well-matched to the mid-IR seeing at the VLT site 
under reasonable optical seeing conditions; thus, VISIR A^-band 
imaging is usually diffraction-limited. A sample of about 30 lo- 
cal AGN were imaged in the two mentioned papers. Due to the 
unprecedented angular resolution, the unresolved core fluxes are 
the best (i.e. least biased by extra-nuclear emission) ones avail- 
able for these objects. In a follow-up work, Horst et al. (2009, 
hereafter Paper III) compared VISIR multi-filter photometric 
spectral energy distributions (SEDs) with sensitive Spitzer IRS 
spectra of the same sources and quantified the level of contami- 
nation by extended emission, finding non-negligible contamina- 
tion in Spitzer observations of about 40% of cases. 

By extracting intrinsic X-ray luminosities (Z/j-io ke\s here- 
after Lx) from the recent literature, we were able to compare 
these with the IR powers measured and found a good correlation 
between these quantities (see, e.g., Fig. 1 in Paper II). The main 
results of our works are that i) the dispersion in the mid-IR:X- 
ray correlation is smaller than found in previous works and ii) 
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both Seyfert-ls and Seyfert-2s closely follow the same correla- 
tion. This second result was not predicted by radiative transfer 
models in which intrinsic AGN emission is re-processed by ob- 
scuring tori with smoothly-distributed dust. This is because an 
optically-thick line-of-sight through the torus will preferentially 
show cooler dust and a lower mid-infrared luminosity for the 
same Lx as compared to an optically-thin one (see, e.g., Fig. 1 
of Pier & Krolik 1993). Result ii) above was also found in pre- 
vious studies with lower-resolution ground and space observato- 
ries (e.g., Krabbe et al. 2001; Lutz et al. 2004), but was ascribed 
to the expected differences being mitigated due to the presence 
of unresolved contamination and other selection biases. On the 
other hand, we found a similar result in our high-resolution 
VLT imaging, even though contamination is much less severe. 
This suggested, instead, that dust clouds have significant dumpi- 
ness in the AGN tori (e.g., Dullemond & van Bemmel 2005; 
Honig et al. 2006; Nenkova et al. 2008a,b; Schartmann et al. 
2008) which leads to the observed mid-IR emission from hot il- 
luminated cloud faces in the inner torus regions being largely 
isotropic. The extent to which the mid-IR emission is isotropic 
can be determined by quantifying the scatter and offsets of ob- 
scured and unobscured AGN at any given Lx in the mid-IR:X- 
ray correlation. Our study was so far restricted to sources with 
low obscuring column densities of cold gas. But the differences 
between the expectations of clumpy and smooth torus distribu- 
tions increase with line-of-sight optical depth, so it is important 
to obtain high-resolution mid-IR imaging of highly obscured 
sources. 

Even irrespective of the question of torus physical structure, 
an accurate census of local Seyferts in the mid-IR is important 
in order to determine their Bolometric accretion energetics and 
to disentangle AGN from starburst emission. If, after removal 
of all selection biases, the scatter in the correlation between IR 
and intrinsic AGN power can be reduced, then the hope is that 
mid-IR observation may be an independent isotropic indicator 
of AGN energetics, similar to forbidden [Om] /15007A emission 
line, far-infrared, hard X-ray or low-frequency radio observa- 
tions (cf. Mulchaey et al. 1994). 

In this paper, we present new high-resolution mid-IR fluxes 
around 12-/im of 16 local Seyferts, including 7 sources ob- 
scured by high Compton-thin absorption (with Vh>10 23 cm -2 ) 
and 6 Compton-thick ones (V H >1.5xl0 24 cm" 2 ). Adding litera- 
ture observations of another well-studied Compton-thick source 
NGC 1068, as well as our previous sample of 25 detections from 
Paper II, and comparing these unbiased core mid-IR measure- 
ments with the best Lx values from the literature, we find the 
luminosity relationship to be very tight, especially for 'well- 
resolved' targets. Various implications for the core structures of 
nearby AGN are discussed and useful broad-band (mid-IR to X- 
ray) spectral relations are derived. A flat Universe with a Hubble 
constant Ho = 71 km s -1 Mpc -1 , Oa = 0.73 and Qm = 0.27 is as- 
sumed. All literature measurements have been converted to this 
Cosmology. 

2. Observations 

A sample of 16 targets was observed with VISIR during ser- 
vice nights spanning the semester from 2007 October to 2008 
March. These were mainly selected from sources that had been 
observed (or had planned forthcoming observations) with hard 
X-ray satellites including Suzaku, INTEGRAL and Swift, so that 
their broad-band X-ray spectra could be modeled and 2-10 keV 
AGN luminosities accurately determined; i.e. the sample is not 
complete. The sky conditions were required to be clear with 0"8 



optical seeing or less, during ESO period 80. Observational and 
data reduction procedures are very similar to those described in 
Paper II; only minimal details are stated here. A standard paral- 
lel chopping and nodding technique with a chop throw of 8" was 
used. Observations were carried out with the small field camera 
(0"75/pixel) in one or more narrow and intermediate-band filters 
covering the TV-band, for each source. As we are interested in de- 
termining source continuum luminosities, filters were selected to 
avoid emission lines and to include continuum regions free of the 
Silicate absorption trough, depending on source redshift. 

The science and standard photometric star calibration frames 
were reduced with a dedicated pipeline (Pantin et al. 2005), in- 
cluding removal of detector features and background subtrac- 
tion. The count rate for one full exposure was calculated as the 
mean of the 3 beams from all nodding cycles of this exposure. 
In order to minimise the effect of residual sky background, we 
chose relatively-small apertures (=s 10 pixels = 1"27) for pho- 
tometry of most sources and corrected the obtained count rates 
using the radial profiles of standard stars. For a handful of faint 
targets, profile-fitting (based on the profiles of bright standard- 
stars) on binned images was carried out instead. As an initial 
error estimate, the standard deviation of the fluxes amongst the 
three nodding beams was used, except for a few cases where the 
statistical error on source counts turned out to be larger. The ab- 
solute flux uncertainty is dependent on the counts-to-flux 'con- 
version' factors determined from standard star observations and 
the known accuracy of stellar templates. For the PAH2 filter, the 
conversion factor has been found to be constant to within « 5% 
over a whole observing period, and this minimum uncertainty is 
added in quadrature to the initial error estimate for all sources. 

The X-ray data used in this paper have been carefully se- 
lected (or computed) from recent published studies including 
spectral modeling over broad energy ranges where available. 
The Lx values are corrected for obscuration, complex spectral 
components and redshift effects, and should be a measure of 
the true, intrinsic 2-10 keV AGN power responsible for heat- 
ing torus clouds. Of course, in cases of extreme obscuration or 
complexity, these corrections can have large uncertainties - see 
§ 6. In this respect, it will be useful to compare against hard X- 
ray luminosities measured above 10 keV (cf. Mushotzky et al. 
2008) in future work, as these should be less susceptible to such 
complexities. Since data from a variety of published sources is 
compiled for each source, the error assignment on Lx is nec- 
essarily somewhat loose. As in Papers I and II, the uncertainty 
is based on the range of peak-to-peak published values, which 
should include source variability and other biases. Each source 
is detailed in the Appendix. 

3. Results 

All 16 sources were detected in filters around 12-yum, and the 
images were dominated in each case by a single IR AGN core. 
The presentation of images and discussion of full multi-filter N- 
band properties is not the purpose of the present work. We re- 
fer any reader interested in the general high spatial-resolution 
Seyfert morphologies to our Paper III, which includes our en- 
tire previous sample. For the purpose of the IR:X-ray luminosity 
correlation (cf. Papers I and II), the 12.3-yum source flux is the 
important quantity. We therefore concentrate mainly on the filter 
most appropriate for probing the continuum closest to 12.3 fan. 
Table 1 lists these filters (referred to as primary filters), the cor- 
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responding source fluxes as well as monochromatic luminosities 
at the central filter wavelength (ALx, hereafter Lmir). We assume 
a flat mid-IR SED in these units, i.e. no ^-correction to 12.3 //m; 
instead, AL 123 of a source which has an SED flat in L\ units, 
say, will be underestimated by only 12 per cent (0.05 dex) in the 
worst case. In addition to our measurements, we included liter- 
ature fluxes of NGC 1068 (Galliano et al. 2005b), as it is one 
of the best-studied of nearby AGN. Some secondary filter fluxes 
are also listed in the table, which we return to in § 5.3. 

Figure 1 shows the correlation between Lx and Lmir- All 25 
VISIR detections from Paper II and the 17 measurements from 
Table 1 are plotted and split into Seyfert 1-1.5 (12 objects, here- 
after 'Sy 1'), Seyfert 1.8-2 or lh (19 objects, hereafter 'Sy 2'), 
LINERs (3 objects) and Compton-thick (8 objects) AGN. The 
' lh' sub-class of Sy 2s refers to sources with broad optical emis- 
sion detected in polarized light (Veron-Cetty & Veron 2006). 
Since Compton-thick sources (with A^h>1-5x10 24 cm -2 ) are of 
special interest, we chose to plot them with a separate symbol, 
though all of them fall under the optical 'Sy 2' category. In ad- 
dition to these 42 sources, one additional Compton-thick object 
(ESO 428-G014) is plotted with a box around it. This is not part 
of the correlation fit, but is used an an aposteriori check; it will 
be discussed in § 6. 

As discussed in Paper II, and shown in Table 1, the sample 
covers a wide range (about a factor 40) of spatial scales at the 
mid-IR angular resolution limit (6q) of the VLT. In order to make 
a physically-meaningful selection, Horst et al. (2008) chose only 
those sources in which the resolved scale (ro), in units of the 
size of the dust sublimation radius (r su b), lay below a certain 
threshold. By comparing Lmir to Lx over the range of resolved 
scales, they identified a value ro/r su b = 560 beyond which the 
ratio Lmir/Lx exhibited a sharp increase in its mean value and 
dispersion. A constant #o=0"35 was assumed, and the relation 

r sub = 0.5^L^pc (1) 

was used, where L 4 ^ is the source Bolometric luminosity in 
units of 10 45 erg s -1 . This relation is approximately correct for 
a sublimation temperature of 1500 K of graphite grains with an 
average radial size of 0.05 yum (e.g., Barvainis 1987, Kishimoto 
et al. 2007; we have also assumed that the bulk of the bolomet- 
ric contribution is due to ultraviolet photons). All sources with 
ro/rsub l ess than this threshold were classified as being 'well- 
resolved.' In other words, the threshold is a simple multiple 
of the expected physical size of the torus, and the suggestion 
is that the rest of the iess-resolved' sources are likely to have 
non-negligible contamination from extended emission that falls 
within the central diffraction Airy ring on the VISIR images. 

In this paper we use the same source classification criterion 
for consistency; we will discuss its limitations and the effects 
of different assumptions in § 5. We thus find 22 (of the total 42) 
sources to be well-resolved, including all Compton-thick sources 
except NGC 5135. Performing a least-square fit on this sam- 
ple with a simple power-law (i.e. a linear relation of the form 
logLviiR= a + MogLx in log-space) gives the following best-fit 
correlation: 

!°g(^§) = (0-19 + 0.05) + (1.11 ±0.07) log(^) (2) 

where the luminosities are normalized to 10 43 erg s _1 , the ap- 
proximate mean power of the sample. The regression accounts 
for luminosity errors on both axes and the uncertainties are based 
on Ax 2 = 1 for each best-fit parameter (cf. fitexy procedure 



in Press et al. 1992). This correlation curve is over-plotted in 
Fig. 1. Eq. 2 is completely consistent with (but slightly tighter 
than) that derived from the smaller sample in Paper II. For easy 
comparison with that paper, we also state the correlation fit with- 
out any normalization, in units of erg s -1 directly: logLMiR=(- 
4.37±3.08)+(1.106±0.071)logL x . The formal significance of 
the correlation, according to the well-known Spearman Rank co- 
efficient (p), is p = 0.93 with a standard Student-f null signifi- 
cance level of 3 x 10~ 10 . As discussed in Paper II, the correlation 
is not just a consequence of correlated distances on both axes; the 
(intrinsic) flux-flux correlation is strong as well, with p = 0.86 at 
a null significance level of 3 x 10~ 7 . This can also be seen from 
a partial correlation analysis between logLx and logLviiR, with 
the correlation component due to redshift removed; these values 
are listed in Table 2. 

When fitting the correlation to all 42 sources without any se- 
lection, the change in parameters is closely equivalent to a sim- 
ple increase in the intercept on the axis of logLMiR by a value 
of 0.24 (i.e., a mean Lmir larger by a factor of 1.7, consistent 
with inclusion of more IR flux in the cores of the less-resolved 
sub-sample); see dotted line in Fig. 1. In this case, p = 0.88 at a 
null significance level of 1 x 10~ 14 . On the other hand, exclud- 
ing the Compton-thick sources has an insignificant effect. A few 
representative fits to various sub-sample populations are listed in 
Table 2. 



4. The intrinsic mid-IR emission of AGN cores 

4.1. Dispersion in Lmir vs. Lx 

Lmir correlates tightly with Lx over three orders of magni- 
tude in Seyfert luminosity. The small scatter is emphasized 
by plotting the distribution of the log-space luminosity ratio 
r=log(LMiR/Lx) in Fig. 2. Our full dataset of 42 sources has a 
mean and standard-deviation of (r, ov)=(0.30, 0.36) in this ratio, 
while the corresponding values for the 22 well-resolved sources 
of all types are (0.15, 0.23). The 15 Compton-thin well-resolved 
sources have the least scatter of all with corresponding values of 
(0.17, 0.18); see Table 2 and right-hand plot in Fig. 1. The disper- 
sion in logLMiR with respect to the best-fit correlation curves (i.e. 
observed vs. predicted logLviiR) in the individual sub-samples 
are very similar to the above 07 values. 

The light blue histogram in the same figure refers to the ISO 
study of Lutz et al. (2004), for which we present their distribu- 
tion in log(/lL^ m /Lx) of 49 AGN - the ones with 6-p.m contin- 
uum detections for which intrinsic X-ray flux estimates could be 
determined. We also removed any sources likely to have signif- 
icant IR starburst contamination (as defined by these authors). 
For this distribution, (r, 07) = (0.53, 0.45). The spread is compa- 
rable to (only slightly larger than) the scatter for our full sample 
of 42 sources. The mean value, on the other hand, is significantly 
larger than for our sample, which may be either due to residual 
contamination within the larger ISO apertures or to the heteroge- 
neous nature of their sample; we return to this in § 4.6. We can at 
least confidently rule out this difference being due to their bluer 6 
p.m central wavelengths. Using median mid-IR AGN SEDs pre- 
sented by Hao et al. (2007), we find that ALf m for Sy 1 (Sy 

2) is about 19 (54) per cent lower than AL l23/im , which is all 
due to red continuum shapes and not due to contamination by 
PAH emission or Silicate absorption. Thus wavelength differ- 
ences alone should have resulted in a lower mean \og(AL^ im /Lx) 
than that for our sample. If we do convert the Lutz et al. lumi- 



4 



P. Gandhi et al.: Resolving the mid-infrared cores of local Seyferts 



45 



^r-. 44 - 
\n 

CD 
CD 



43 - 



> 

CD 
O 

I 



^ 42 P 

CD 
O 



41 - 



■ Sy 1-1.5 
♦ Sy 1.8-2 
C Compton-thick 
▲ LINERs 




41 42 43 44 

\ogXL x (12.3 |im) [erg s" 1 ] 



45 



45 



„ 44 

to 

CD 
1— 

? 43 



^ 42 1 1 

CD 
O 



41 - 



■ Sy 1-1.5 

♦ Sy 1.8-2 
▲ LINERs 




41 42 43 44 

\oqXL x (12.3 |xm) [erg s" 1 ] 



45 



Fig. 1. (Top): Mid-infrared:X-ray luminosity correlation for the sample of 42 VISIR detections from Table 1 and Paper II. One additional object 
plotted (to make 43 in all) is ESO 428-G014 (the 'C with both a box and a circle around it), which is not part of our fit sample, but used to validate 
the correlation in the Compton-thick regime (see § 6). The sources classified as being 'well-resolved' are encircled, and the diagonal dashed line 
is the fitted correlation for these (Eq. 2 and Table 2). The dotted grey line is the correlation fit to all 42 sources. Bottom): Lmir-vs.-Lx for the 34 
Compton-thin sources only, with the correlation fitted to 15 well-resolved sources (Table 2). 
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Table 1. Observational parameters of our sample of 16 VISIR detections plus archival data on NGC 1068, tabulated in order of increasing Right 
Ascension. (2) Seyfert type, mostly according to Veron-Cetty & Veron (2006), except for the two Swift targets, ESO 323-G032 and NGC 4992; 
see Appendix. (3) Redshift z. Values marked d are instead distances in Mpc. (4) Obscuring gas column density in units of log [cm -2 ] taken from 
literature X-ray spectral fits. (5) Approximate resolved scale limit in pc, for an angular resolution of 0"35. (6) Name and central wavelength in 
/jm of a secondary filter chosen to test systematic variations in the L M i R -L x correlation (details in § 5.3). (7) Secondary filter flux. (8) Chosen 
primary filter name and central wavelength in /im. *For NGC 1068, the fluxes are from the high angular resolution (0"2-0"4) Keck and Subaru 
studies (Mason et al. 2006; Tomono et al. 2001) at slightly-differing wavelengths. (9) Observed source core flux in primary filter. (10) log of 
monochromatic primary filter luminosity expressed in erg s~' . For our mid-IR:X-ray correlation we assume a flat SED in AL^, i.e. no ^-correction 
to 12.3 /im. (11) log of 2-10 keV intrinsic X-ray luminosity. See Appendix for details on the literature used and notes on individual sources. 



Sub-sample 


N 


P-z 


a 


b 


r 


o> 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


All 


42 


0.81 


0.41+0.03 


1.11+0.04 


0.30 


0.36 


All, excluding CT 


34 


0.79 


0.43+0.04 


1.11+0.04 


0.35 


0.36 


All Sy 1 


12 


0.71 


0.34+0.09 


1.18+0.14 


0.34 


0.31 


All Sy 2 + CT 


27 


0.77 


0.40+0.06 


1.40+0.11 


0.25 


0.38 


Well-resolved 


22 


0.89 


0.19+0.05 


1.11+0.07 


0.15 


0.23 


Well-resolved, excluding CT 


15 


0.93 


0.18+0.05 


1.08+0.07 


0.17 


0.18 


Less-resolved 


20 


0.76 


0.62+0.05 


1.11+0.06 


0.47 


0.40 




Systematic variations 








Well-resolved (Systematic IT) 


25 


0.82 


0.25+0.05 


1.21+0.08 


0.17 


0.30 


Well-resolved (Systematic 1J,) 


16 


0.90 


0.26+0.05 


1.07+0.08 


0.21 


0.23 


Well-resolved (Systematic 2) 


22 


0.83 


0.22+0.05 


1.09+0.07 


0.16 


0.27 


Well-resolved (Systematic 3) ' 


* 


* 


1.87+0.31 


1.11+0.10 


* 


* 


Well-resolved (Systematic 4) 


* 


* 


0.15+0.05 


1.02+0.07 


* 


* 


All (Systematic 5) 


42 


0.74 


0.42+0.03 


1.11+0.04 


0.30 


0.44 


Well-resolved (Systematic 5) 


22 


0.78 


0.21+0.05 


1.09+0.07 


0.13 


0.31 



Table 2. Correlation properties between logLx and logLMiR to various sub-sample populations (Column 1), each containing N sources (Col. 2). 
The partial correlation coefficient (p. z ), with the effects of redshift removed are listed in Col. 3. Regression of the form logL^ IR = a+blogL^ 

(luminosities in units of 10 43 erg s~') was carried out; the best-fit parameters are listed in Cols. 4 and 5 along with their 1-tr uncertainties. Cols. 6 
and 7 list the mean value and scatter in r=log(LMiR/ix); see § 4.1. 'CT' refers to Compton-thick sources. The two solutions likely to be of most 
interest (see § 7) are written in bold. The final seven rows refer to tests for systematic variations in the correlation; see § 5.1- 5.3. The test fit 
marked with a ' (Systematic 3) is carried out in linear (not log) luminosity space (see § 5.2). For Systematics 3 and 4, the sub-samples are identical 
to the main well-resolved sample of 22 sources, and all '*' parameters are identical as well (see bold text under 'Well-resolved'). 



nosities to rest-frame AL A ' M using the Sy 1, Sy 2 and quasar 
fe-corrections from Hao et al., we find (r, 07) = (0.71, 0.50). 

With regard to the smaller dispersion in our well-resolved 
sub-sample, is this simply due a serendipitous selection of a few 
sources from our full sample which, in reality, have no physi- 
cal inter-relation? The dispersion of a smaller sample may well 
be measured to be smaller due to statistical, rather than phys- 
ical reasons, as pointed out by Lutz et al. (2004). We rule out 



this possibility based on two arguments. First, note the clear 
and large differences between the scatter of the well-resolved 
sub-sample and of other sub-samples containing approximately 
equal objects in Table 2. The difference is especially clear when 
comparing against the less-resolved sources which have much 
larger (f, cr r )=(0.5, 0.4) values. Next, in order to further esti- 
mate the effect of any fortuitous selection, we constructed 30000 
bootstrap populations of 22 sources each, drawn randomly from 
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the full sample of 42 sources without replacement. We measured 
the dispersion in each of these populations and computed the 
fraction of random populations with a dispersion less than or 
equal to the dispersion measured in the data. The result is that 
our measured data dispersion is robust to a serendipitous selec- 
tion at 99.2 per cent; only about 250 random sub-samples had a 
variance smaller than that of the data. 

4.2. The compositions of the unresolved mid-IR cores 

The imaged VISIR cores of the vast majority of our full sample 
appear completely unresolved, even at the high-resolution of the 
VLT (as mentioned above; see also extensive discussion in Paper 
III). The physical resolution limits (ro) of our observations range 
over tens to hundreds of pc (see Table 1; Paper II). This is larger 
by factors of several hundred than an inner torus radius com- 
puted according to Eq. 1 . The full radial extent of the torus, on 
the other hand, is more difficult to constrain (e.g. Granato et al. 
1997). But our resolution limit is still larger by factors of a few 
to several tens than the outer radii suggested by recent clumpy 
torus models (cf. Honig et al. 2006, Nenkova et al. 2008b and 
references therein). So our listed ro values should be taken as 
strong upper limits of the physical sizes of the individual tori (at 
least the size relevant for the distribution of hot mid-IR-emitting 
dust). 

What about the composition of the unresolved AGN cores? 
The small scatter in our correlation can be used to constrain their 
average composition to some extent because any contribution by 
multiple, unresolved components is bound to increase the dis- 
persion in log(LiR/Lx)- 

A very tight correlation may argue for a single population 
of particles emitting both X-ray and IR photons via some non- 
thermal process. Though we cannot discount an important non- 
thermal mid-IR component for every one of our sources, such a 
component is not the dominant cause of the tight correlation, 
as can be argued in several ways. Only a few of our targets 
show significant non-thermal emission in the radio, including 
3C 445 and Cen A (in which a non-thermal mid-IR synchrotron 
component has indeed been inferred to exist by Meisenheimer 
et al. 2007; though see Radomski et al. 2008 for other possi- 
bilities). Furthermore, the X-ray power-laws, if extended to the 
IR regime, severely underestimate the observed VISIR fluxes 
in most cases (but see Carleton et al. 1987). Finally, extensive 
modeling of individual sources also supports the weakness of 
any non-thermal component. An extended mid-IR structure has 
been resolved out with interferometric observations (Jaffe et al. 
2004) in at least one of the objects, NGC 1068. Honig et al. 
(2008a) have modeled the broad-band emission from this source 
and placed strict limits on the non-thermal radiation, favouring 
instead thermal emission from a clumpy torus cloud medium. 

In any case, our empirical luminosity correlation is valid and 
important, regardless of its underlying physical origin. We leave 
further discussion of this issue to future works including de- 
tailed monitoring observations and modeling of AGN tori; in- 
stead, we focus now on another consequence of the small cor- 
relation spread - limits on the spatial extent and power of star- 
formation (SF) in the mid-IR at high spatial resolution. 

4.3. Constraining unresolved mid-IR star-formation 

Powerful nuclear starbursts (SB) should produce copious 
amounts of infrared emission, but comparatively little X-rays 
(e.g. Rodriguez Espinosa et al. 1987; Ranalli et al. 2003). While 



the bulk of this is radiated at far-IR wavelengths (at ^60-yum), 
the mid-IR flux around 10-yum due to SF is far from negligible 
(e.g. Maiolino et al. 1995, Galliano et al. 2005a, Siebenmorgen 
et al. 2008, our Papers II and III). Furthermore, the instantaneous 
IR luminosity of a SB can vary by more than an order of mag- 
nitude (e.g. Cid Fernandes et al. 2004; Davies et al. 2007), de- 
pending on its size and age. The resultant expected scatter in the 
combined AGN+SB power can be large. On the other hand, our 
tight mid-IR:X-ray correlation suggests that AGN emission in 
the point cores of our high-resolution imaging completely dom- 
inates over any SF on similar scales at mid-IR wavelengths, and 
provides us the opportunity to place limits on any unresolved SB 
contribution. 

In a recent work, Ballantyne (2008) investigated the efficacy 
of nuclear SBs at z ~ 1 to obscure the sources that power the 
bulk of the cosmic X-ray background radiation. Although the 
physical motive and cosmic epoch of Ballantyne's study is dif- 
ferent from ours, his numerical models make useful predictions 
of the emitted mid-IR (from the SB) and X-ray (AGN) fluxes, 
which can be compared with our observations. In particular, the 

expected distribution of log(/LL^ m /Lx) from his Fig. 1 1 is il- 
lustrated with the green histogram in Fig. 2. For this modeled 
distribution, cr r a:0.75, a factor of about 3 (2) greater than our 
sample of well-resolved (all) sources, though it must be noted 
that the nuclear SB luminosity ratios are indicative values only, 
according to Ballantyne. Conversion of the model luminosities 
to 12 fim should not affect the dispersion much (cf. Fig. 7 in 
Ballantyne 2008, in which the 6 and 12 /im SB luminosities are 
approximately equal). 

SF is certainly present around the cores of Seyferts, but it is 
the spatial scale being probed that is the important quantity here. 
As an illustrative example of this, in Fig. 3 we show our VISIR 
high-resolution imaging of one target - IC 3639 - whose nuclear 
SB was studied in detail by Gonzalez Delgado et al. (1998) and 
found to have a bolometric luminosity (determined mainly from 
ultraviolet observations) very similar to the total power emitted 
by the central AGN. The point-like mid-IR profile is in sharp 
contrast to the over-plotted contours of HST imaging obtained 
by these authors. Images of all our targets from Papers I and II 
can be found in Paper III (Horst et al. 2009); other multi-band 
images and SEDs from the present sample will be presented in a 
future work. 

The median resolved scale (ro) for our 22 targets classified as 
being well-resolved is 70+ 1 1 pc, and for the rest of the 20 less- 
resolved targets is 1 15 ±24 pc (errors correspond to the standard 
deviation of the median of 10000 bootstrap samples drawn with 
replacement). Another way to gauge the difference between the 
two samples is a two-sided Kolmogorov-Smirnov (KS) test, ac- 
cording to which the distributions of resolved scales differ at the 
97 per cent level. On the other hand, the distributions of r su b val- 
ues of the two samples are consistent with each other. In other 
words, the dominant factor governing source classification turns 
out to be simply the distance, with the well-resolved sources be- 
ing on average closer to us (their median redshift is 0.010, while 
Zmedian for the rest is 0.017) and with smaller ro values. 

It then follows that the relevant scale within which the AGN 
completely dominates the SB is the central ~ 60 - 80 pc (the me- 
dian of the well-resolved sample) of the nuclei of local Seyferts. 
In Table 2, we found o> = 0.23 for the well-resolved sample. 
If this scatter were to be solely due to nuclear SB and assuming 
negligible X-rays from them, it would imply that SB within these 
physical scales typically have a 12 /im luminosity of no more 
than l-10 -ov = 40% of the total unresolved (AGN+SB) L M1R . 
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Fig. 2. Dispersion in logLMiR/^x- For our data, three distributions are plotted in order of increasing depth into the figure: i) that for 15 well-resolved 
Compton-thin sources is shown by the pink cross-hatched histogram; ii) that for all 22 well-resolved sources by the black hatched area; and Hi) 
that for all 42 sources by the unfilled histogram with dashed outline. The grey filled circle and arrow shows the mean and extreme values of the 
correlation curve in Eq. 2 over the fitted luminosity range. For the sample of Lutz et al. (2004), the light-blue hatched area shows the distribution 
of 49 sources in log (AL 6 ^' m /Lx), normalized to the maximum value of our full data histogram. Finally, the light-green, dotted histogram is the 
model prediction for the AGN+starburst model of Ballantyne (2008), also normalized. 



Assuming the dispersion for the full sample of 42 sources in- 
creases this limit to 56% while the use of the Compton-thin sub- 
sample reduces it to 34%. These are conservative upper-limits, 
because other effects such as AGN infrared variability and the 
scatter due to statistically-differing distributions of torus clouds 
will introduce additional uncertainty; so the true constraint will 
undoubtedly be even tighter. 

Beyond a radius of ~ 100 pc (using the less-resolved sam- 
ple alone), the limit on the SB power increases to ~ 1 .5 times 
the 12 yum luminosity of the AGN itself, allowing for the bulk 
of star-formation to occur on these scales. This is completely 
consistent with other recent studies: e.g. Imanishi (2002), Soifer 
et al. (2003) and Davies et al. (2007). The last authors ob- 
tained detailed integral-field near-IR spectroscopy of several lo- 
cal Seyferts and found AGN and nuclear SB powers to be com- 
parable only on kpc scales, /f-corrections of their results to 12- 
yum are expected to be small compared to the differences due to 
integration over different spatial scales. Finally, we emphasize 
that although our observed scatter is much smaller than that pre- 



sented by Ballantyne (2008) this is not, by any means, a negation 
of their model which is relevant at z ~ 1 . It simply means that we 
are resolving out any similar SBs at the present epoch, if they 
exist. 

In summary, we suggest that our high-resolution observa- 
tions are predominantly probing the nuclear torus emission, un- 
biased by any overwhelming contamination due to powerful SF. 

4.4. lnfrared:X-ray spectral indices 

Having determined the true mid-IR and X-ray luminosities at- 
tributable to the AGN themselves, we can now derive their in- 
trinsic SED shapes. A spectral index connecting two widely- 
separated frequencies is often used to describe sparsely sam- 
pled broad-band spectra. In accordance with the widely-used 
optical:X-ray spectral index (aox), we define the index 



2 keV'^12.3 urn' 
log(v 2 kev/vi2.3 fim) 



: 0.233 log(^ keV /L* ) 



(3) 
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Fig. 3. VISIR NEII.2 (Ac = 13.04-yian) image of IC 3639, showing the 
clean point-like core profile. North is up and East is to the left; the 
image was obtained at 0"75/pixel sampling and is 2" x 2" wide. The 
overlaid contours are from the ultraviolet (UV) Faint Object Camera 
image obtained from the HST archive. Contour levels range from 10% 
to 100% of the peak UV flux. The extended SB emission over the central 
arcsec, especially the central 'horizontal-S' shaped structure of knots 
(cf. Gonzalez Delgado et al. 1998) are clearly visible in the UV, but not 
at 12-yum. 



where U represent the monochromatic luminosities [erg s 
Hz -1 ] at frequencies (v) corresponding to 2 keV and 12.3 /im, 
both in the rest-frame. We have computed IZ keV from the indi- 
vidual X-ray photon-indices (usually determined from spectral 
fits over the ~ 0.5 - 12 keV range) in the literature. 

We then find a mean arx=-l .10+0.01 (0.06) and -1 . 1 8+0.02 
(0.10) for the well-resolved, and less-resolved, samples respec- 
tively, and -1.14+0.01 (0.08) for the whole sample; the num- 
bers in brackets list the full scatter (cr,^). This distribution is 
shown in Fig. 4. The plot shows that not only is the spectral slope 
for the less-resolved sources steeper (consistent with including 
more mid-IR emission for the same Lx), but that the dispersion 
of slopes for these sources is larger. 

Often in deep surveys, it is not possible to determine the X- 
ray photon-index from individual spectral fits, due to the faint- 
ness of sources. For such cases, we note it is a good approxi- 
mation to assuming a constant r w 1.9 over the 2-10 keV range. 
This is because the typical observed range of F is » 1 .5-2.3 (e.g. 
Nandra & Pounds 1994; Mateos et al. 2005; Piconcelli et al. 
2005), and the difference between a constant r = 1 .9 and the 
extremes of the stated range results in only a small difference 
of Ao'ix=+0.04 (assuming no change in Lx). In terms of pre- 
dicting the luminosity of a source at one frequency from a flux 
(i.e. luminosity) measurement at the other frequency, using the 
average spectral slope computed above will give an answer accu- 
rate to within a factor of 1.5, over the mentioned range of X-ray 
photon-indices. 

It is well-known that the ultraviolet:X-ray spectral index 
(aox) of AGN is strongly dependent upon luminosity (e.g. 
Vignali et al. 2003; Steffen et al. 2006). The correlation of Eq. 2 
is only slightly non-linear, which implies that the correspond- 
ing mid-IR:X-ray index should instead depend only weakly on 



luminosity. Assuming again that a constant X-ray photon-index 
(F = 1 .9) is applicable to all sources for simplicity, the expected 
power-law index of dependence of a\x on the monochromatic 
IR luminosity is easily derived from Eq. 2 and Eq. 3 to be 
0.233(£> -1 - 1), where b is the best-fit correlation slope. We found 
b = +1.11 for the well-resolved sub-sample, which implies an 
expected dependence of a\x °c (L^j IR )~ ' 02 . The right-hand panel 
of Fig. 4 shows the full dataset of individual sources in the crrx- 
vs-Lj^ IR plane; the errors on arrx incorporate the uncertainties on 
both X-ray and mid-IR luminosities. The dotted line shows the 
regression of the well-resolved sources, and its best-fit slope of 
-0.022+0.013 matches the expected value. 

4.5. Bolometric corrections 

Correction factors (ffloi) for obtaining full bolometric AGN 
power emissions based on observations over limited energy 
ranges are crucial for computing the total accretion energy re- 
lease in the Cosmos, as well as studying black hole growth 
and evolution (e.g. Elvis et al. 1994; Fabian & Iwasawa 1999; 
Marconi et al. 2004; Hopkins et al. 2007; Gilli et al. 2007). In 
particular, it is natural to connect source distributions seen at 
mid-IR and X-ray energies (e.g. Gandhi & Fabian 2003) as these 
are also good probes of any evolution in obscured AGN popula- 
tions (e.g. Ballantyne & Papovich 2007; Treister et al. 2008). 

Intrinsic bolometric corrections connecting the absorption- 
corrected X-ray power with the integrated (optical/ultraviolet to 
hard X-ray) emission have been found to increase with AGN lu- 
minosity (e.g. Marconi et al. 2004; see also above discussion on 
ffox)- If the mid-infrared flux, on the other hand, is the result of 
reprocessing and thermalization of intrinsic AGN photons within 
torus dust clouds, then the infrared correction factor will depend 
upon radiative transfer within each cloud, as well as the clouds 
covering factor. 

From an empirical standpoint, the expected infrared bolo- 
metric corrections (fg^^Boi/LMiR) can be easily derived by 
combining our Lmir : Lx correlation with the Lq \ : L2-10 trend 
found by Marconi et al. (2004). These last authors present the 
X-ray bolometric corrections (fg ol =LBoi/L2-io) as a function of 
the scaled bolometric luminosity X.-log(Lg \/LQ)-12, and we 
adopt the same units to ease comparison. The value of Lq used 
is 3.826xl0 33 erg s -1 . The following equation is obtained: 

!°g f Bof (-0 = W°g f Boi(^) -« + d- b)(X + 12 + logL - 43) (4) 

Expressing it in this form makes clear the dependence on our 
fitted slope and intercept (which can be picked from Table 2 for 
any of the fits), and the additional log-luminosity units collected 
at the end. It also aids in Monte Carlo determination of confi- 
dence intervals, because Marconi et al. (2004) state an uncer- 
tainty of 0. 1 in logfg ol and we have computed the uncertainties 
in a and b in Table 2. Assuming these to be normally-distributed, 
the resultant 1-cr uncertainty can be derived from the scatter in 
a large number (we used 10000) of randomizations of Eq. 4 at 
various values of Substituting the best-fit a and b for our well- 
resolved sample (Eq. 2), and expanding logfg ol from the relevant 
Eq. 21 in Marconi et al. (2004) then gives 

lo g f Boi R (^) = l-24 + 0.16i: + 0.012i: 2 -0.0015X 3 (5) 

This curve is drawn in Fig. 5, along with the computed uncer- 
tainty. The plot shows a monotonically-increasing fg 4 ^ correc- 
tion for the Sy luminosities that we probe, with a range of av- 
erage values of « 10 - 30, while the X-ray correction is steeper. 
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Fig. 4. (Left): The two-point broad-band SEDs of the well-resolved (purple, solid lines) and less-resolved (green, dashed lines) sub-samples, 
showing the distribution of ayx values. Within both sub-samples, sources have been normalized for clarity to the average sub-sample luminosity 
at the mean frequency of logv=15.54. The thick central lines denote the mean a\x slopes. (Right): The distribution of individual a?ix values as a 
function of log L v l2 3 ^ m . The dotted line is the best-fit regression curve for the well-resolved sub-sample and has a slope of -0.02. 



These values lies at (and above) the upper range of Bolometric 
corrections found by Spinoglio et al. (1995); this is consistent 
with us probing closer in to the individual nuclei than the IRAS 
data used by these authors, but specific differences (e.g. in com- 
putation of LboI, as well as their much larger sample size) need 
to be kept in mind. 

It is important to note a caveat implied by more recent stud- 
ies. Vasudevan & Fabian (2007, 2008) find only a marginal trend 
of fed with Lx, inferring instead a stronger dependence on the 
the Eddington ratio which is proportional to LboI normalized to 
the black hole mass. In this case, infrared bolometric corrections 
may be expected to vary more weakly with luminosity than in 
Fig. 5, or to show more complex behaviour. 

4.6. From Seyferts to quasars? 

Our study has probed local AGN, which are predominantly 
Seyferts. AGN with quasar-like X-ray luminosities (Lx» 10 44 
erg s -1 ) are absent from our sample due to this selection of 
nearby targets and because distant sources are less-resolved. On 
the other hand, it becomes increasing difficult for SF to domi- 
nate the energetics of sources at the high end of bolometric lu- 
minosities. The bulk of the observed IR fluxes of the most lu- 
minous Seyferts and quasars are likely to be due to accretion 
onto the central black hole, irrespective of instrumental angular- 
resolution. Thus, we should be able to check our correlation in 
the luminous quasar regime by using space-telescopes observa- 
tions of distant quasars. 

We have tested this on the sample of Lutz et al. (2004), for 
which the authors have carefully decomposed AGN and host 
galaxy contributions using ISO spectra. The sample includes 
Seyferts, quasars as well as ultra-luminous infrared galaxies 
(ULIRGs), and extends out to z~ 1 . We converted the 6 fim con- 
tinuum luminosities of the 49 sources used in Fig. 2 to 12.3 fim 



using the ^-corrections of Hao et al. (2007) for Sy Is and Sy 
2s below Lx=10 44 erg s -1 , and using the corresponding quasar 
SEDs above this limit; but we note that the conclusions below do 
not depend upon these corrections. The result is shown in Fig. 6. 

It is seen that the ISO sample deviates significantly from our 
best-fit correlation at low luminosities but not at high luminosi- 
ties, just as might be expected due to contamination at low pow- 
ers. The match is particularly good for the brightest unobscured 
quasars (with logLx~44) whose observed emission should be 
the least biased by any extended SF. 

The deviations of the type 2 AGN, on the other hand, are 
large even for luminous objects. The sense of the deviation is 
such that a higher IR power is observed at any given Lx, com- 
pared to the distribution of type 1 AGN. Examining the three 
most luminous type 2 AGN, we find that all of them are actually 
ULIRGs: IRAS 23060+0505 (Brandt et al. 1997), IRAS 05189- 
2524 (Severgnini et al. 2001) and Mrk 463 (Bianchi et al. 
2008). The observed core structure of ULIRGs is quite different 
from Seyferts and quasars. Many ULIRGs are seen to occur in 
strongly interacting and merging systems where gas gets chan- 
neled to the cores and drives very compact and powerful SBs. 
Two X-ray emitting AGN are sometimes seen in close proximity 
- one from each of the merging systems. Therefore, no apriori 
match of ULIRGs (i.e., without extensive corrections) with our 
Sy mid-IR:X-ray correlation is expected. Furthermore, the X-ray 
luminosity of at least one of these - Mrk 463 - has been revised 
upwards by new Chandra imaging (see Bianchi et al. 2008, who 
find a binary AGN with higher Lx than previous estimates), thus 
actually moving it closer to our correlation line. 

With all these caveats in mind, we regard the match of the lu- 
minous sources from Lutz et al. (2004, at least the unobscured 
ones) to be fair; consequently, we suggest that Eq. 2 can be ex- 
tended into the quasar regime as well. Interestingly, our corre- 
lation is well-straddled by the two correlation slopes quoted by 
Fiore et al. (2008, see their Eq. 1) for Sy and quasars respec- 



10 



P. Gandhi et al.: Resolving the mid-infrared cores of local Seyferts 
^12.3 um ( ergs -1) 



10 



42 



1:0 



42 



70 
50 



10 



43 



10 



43 



1Q 44 1Q 45 

1 44 1 45 



L x (erg s 1 ) 



OkeV 




10 



TO" 



10 



TT 



10 



T2~ 



10 



TT 



■^Bol C-Sun) 



Fig. 5. Mid-infrared bolometric correction f^j R =LBol/^MlR ( re d continuous curve) - determined from our Lmir^x correlation for well-resolved 
sources (Eq. 2) and the fg ol =Lg i/Lx trend observed by Marconi et al. (2004). This latter X-ray trend is shown as the blue dashed curve. The red 
hatched area shows the 16th and 84th percentile confidence region in 10000 Monte Carlo randomizations, incorporating normal uncertainties on 
our correlation parameters as well as on logfg o) . 



tively, in the CDF-S and C-COSMOS fields. Of course, further 
detailed validation with enlarged samples will be important. 



5. Systematic effects 

5.1. To be or not to be well-resolved 

As described in § 3, our threshold for classification of sources 
as being well-resolved, or not, is based on locating a discontinu- 
ity in the LMiR/ix-vs.-ro/r su b (luminosity ratio vs. normalized 
scale) parameter space. The reasoning is that for sources at the 
higher end of the distribution of VISIR core physical size scales 
(as compared to the torus), significant unresolved contamination 
is allowed and can lead to an increase in the scatter as well as 
the mean Lmir at any given X-ray power; these sources are the 
'less-resolved' ones; the rest being 'well-resolved'. Due to the 
empirical nature of our search, though, it has important limita- 
tions which we discuss in this section. But we emphasize that the 
correlation for all 42 sources remains robust (and with a small 
scatter relative to previous studies) even if no selection based on 
resolved scale is applied (see Fig. 2). 

Firstly, adding the targets from this paper to the luminosity- 
ratio vs. normalized scale plot of Paper II (Fig. 3 of that work), 
the full new distribution is plotted in the left-hand part of 
Fig. 7. A crude uncertainty in the normalized scale threshold 
can be estimated simply from the mean errors of the binned 
source distribution around the adopted threshold of ro/r su b=560, 
which suggests A(ro/r su b) ~ 40. Thus increasing (decreasing) 
our threshold to ro/r su b=600 (520), the number of well-resolved 
sources changes from 22 to 25 (16). The corresponding best- 
fit slope steepens (flattens) to 1.21 (1.07). The details of these 
fits are listed in Table 2 as 'Well-resolved (Systematic If)' and 
'(Systematic IX)' respectively. 
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Fig. 6. Our correlation from Fig. 1 , compared with the distribution of 
sources from Lutz et al. (2004, cf. their Fig. 5). Their 6 /jm continuum 
luminosities have been converted to 12.3 /im using median ^-corrections 
from Hao et al. (2007). Luminous unobscured AGN (1-1.5; blue empty 
squares; those least biased by contamination) fall smoothly onto our 
correlation in the quasar regime. 



Secondly, in Paper II, a linear scaling of Lbo1=10Lx was as- 
sumed in the computation of r su b. On the other hand, Bolometric 
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corrections may increase with source luminosity (as discussed 
in the previous section), in which case r su b will increase faster 
than implied by the linear dependence above, and will result in 
more of the luminous sources satisfying any given threshold in 
7"oA" S ub- Using the ffioi-^x relation of Marconi et al. (2004), we 
have recomputed r su b (and hence the normalized scale ro/r su b) as 
a function of log(LMiR/Lx); the resultant distribution is shown as 
the light grey points in the right-hand panel of Fig. 7. On aver- 
age, the normalized scales are shifted to lower values due to the 
faster increase in r su b. In order to best determine the scale where 
the average distribution of log(LMiR/£-x) undergoes a change, 
we adopted a floating threshold in normalized scale. The distri- 
butions above and below this floating threshold were compared 
with a two-sided KS test. A value of ro/r su b=420 proved to be 
the best discriminator, with the two distributions (well-resolved 
vs. less-resolved) differing at 98.6 per cent. This can also be 
roughly discerned by eye in the plot as the noticeable change 
in the binned distribution; the mean and scatter of the resultant 
well-resolved and less-resolved sub-samples are (r, cr r )=(0.16, 
0.27) and (0.45, 0.38) respectively. The number of sources which 
satisfy this well-resolved threshold is 22 (a slightly different set 
of 22 from that discussed in all previous sections) and the re- 
sultant Lmir-£x correlation fit is largely unchanged ; see details 
in Table 2 under heading 'Well-resolved (Systematic 2)'. So the 
correlation parameters are insensitive to changes in the resolu- 
tion threshold, as long as the corresponding well-resolved sub- 
sample is self-consistently defined. 

Lastly, we note that while recent work has confirmed the de- 
pendence of r su b°c VZ used in Eq. 1 (Suganuma et al. 2006), the 
absolute normalizing scale (in pc) of the relationship remains 
uncertain. Kishimoto et al. (2007) have shown that the inner- 
most torus radii inferred from near-IR reverberation mapping of 
type 1 AGN are systematically smaller by a factor of ~3 than 
would be implied by the widely-used canonical assumptions im- 
plicit in Eq. 1 . They discuss various alternatives to account for 
this, including a sublimation temperature higher than that used 
for graphite (~1500 K), or a typical grain size of ~2 fim - about 
four times larger than usually quoted. In any case, a smaller r su b 
will result in a simple translation of all our sources to larger x- 
axis (ro/r su b) values in Fig. 7. The large change in log(LMiR/£x) 
will then be seen to occur around 1700 x ro/r su b (using the factor 
~3 offset found by Kishimoto et al.). Again, if the threshold is 
consistently chosen, there is no effect on our correlation. 



5.2. Treatment of errors 

The individual X-ray error bars are a significant source of uncer- 
tainty for the correlation. Our strategy has been to preferentially 
choose published analyses with broad-band spectral modeling, 
and to assign an error based on the range of source variability 
and fit parameters inferred in the recent literature. If the IR is 
indeed due to reprocessing in the torus, then the relevant quanti- 
ties to compare are not necessarily X-ray observations obtained 
simultaneously with the IR, but rather the long-term (few years; 
relevant for a pc-scale torus) X-ray average flux. The IR flux, 
on the other hand, is expected to be much less variable. The 
range of published fits spanning various missions and several 
years (in addition to spectral complexities) results in X-ray error 
bars much larger than the corresponding IR ones. 

For simplicity, we adopted equal-length (high and low) er- 
ror bars in log-luminosity space, which should be approximately 
correct given the many sources of uncertainty. On the other hand, 
it may be more natural to use equal-length error bars in linear 
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Fig. 7. Left (a): Logarithmic luminosity-ratio vs. normalized resolved 
scale (rn/r S ub)> wm ^ r sub being determined according to a constant linear 
scaling of Lg i-vs.-Lx- The primary adopted threshold of ro/r su ^=560 
is shown as the dotted vertical line. Sources falling on the left of this 
line are taken to be 'well-resolved', and the rest 'less-resolved'. All 42 
sources with symbols as in Fig. 1 are shown in grey, with the errors on 
scale accounting for uncertainty in Lx and slight variations in mid-IR 
seeing. Bins in rn/r su b are shown by the dark-green bold circles with 
horizontal error-bars denoting full bin-widths. The standard deviation 
of the mean log-luminosity ratio in each bin is shown by the vertical 
error-bars on these circles. Right (b): In this plot, Lb i (and hence r su b) 
is determined according to the bolometric corrections of Marconi et al. 
(2004). The threshold is the best-fit discriminatory scale at rn/r su b=420. 
See §5.1. 



space because flux (and count) measurement errors are usually 
normally distributed. Thus directly fitting a power-law of the 
form Lu\K-al\ to the base sample of 22 well-resolved sources 
and accounting for the equal-length X-ray errors in luminosity, 
still results in a slope b — 1.11 ± 0.10, i.e. very similar to that of 
Eq. 2; see Table 2 'Systematic 3'. 

Finally, one may question the veracity of the statistical na- 
ture of our errors and ask what happens if we ignore these com- 
pletely. We have tried a fit according to the algorithm of Isobe 
et al. (1990) which can be used when the error assignment and 
the nature of scatter about a linear relation are ill-understood. 
We directly used the log-luminosities from the last two columns 
of Table 1 (and corresponding ones from Paper II) without the 
tabulated errors. Adopting the Ordinary Least-square Bisector 
estimate recommended by Isobe et al. results in a regression 
slope of b — 1 .02 + 0.07, i.e. flatter than previous values and con- 
sistent with unity at l-cr confidence. This procedure effectively 
gives identical weighting to all sources, including highly ob- 
scured ones; but it is precisely these last sources that tend to have 
larger inherent uncertainties (e.g. when correcting for Compton- 
thick obscuration) and that tend to lie above the correlation line 
in Fig. 1 (see discussion that follows in § 5.4 and § 6). We thus 
caution that such a fit may not be wholly appropriate, though 
we include it for completeness. The fitted slope and intercept are 
listed in Table 2 under 'Systematic 4'. 

5.3. Filter selection 

Depending on the redshift of each source, we have tried to 
choose narrow filters with central wavelengths that ought to be 
dominated by continuum emission. But we cannot rule out con- 
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tamination due to emission line wings in every case. For in- 
stance, using the median Sy 2 SED from Hao et al. (2007), we 
estimate a worst-case contamination of the [Nen] il2.81 emis- 
sion line in our NEII_2 filter photometry of about 14 per cent for 
the observation of NGC 3281. 

As mentioned, we observed most sources in a number of fil- 
ters around 12-pm. We can, therefore, gauge any overall bias 
on the correlation fit due to selection of our particular filters by 
choosing instead the second most-appropriate filter from our set, 
i.e. the one whose central wavelength is closest to 12.3-yum, fol- 
lowing the primary filter. These filters, and the related fluxes for 
the new targets in this paper, are listed in Table 1 as 'secondary' 
filters. Including data from Papers I and II, we had access to 
34 secondary filters fluxes out of the 42 sources; for the rest, 
we retain the primary fluxes. Converting these fluxes to 12.3- 
yum luminosities as before, and fitting the new Lmir-^x corre- 
lation reassuringly results in a negligible change of the correla- 
tion fit parameters. Details are listed in Table 2 under the rows 
headed 'Systematic 5'. The dispersion in the correlation (o>), 
on the other hand, does exhibit an increase. This is probably 
a consequence of the extra contamination in some of the bluer 
and broader secondary filters (i.e. SIV, SIC, PAH1; with central 
wavelengths around ~9— 12 yum) due to Silicate absorption, PAH 
and other features. 

In summary, the considerations of the three preceding sub- 
sections suggest systematic variations in the mid-IR:X-ray cor- 
relation slope over ~ 1.02 - 1.21, in addition to statistical fit un- 
certainties. 

5.4. Infrared extinction corrections? 

No extinction corrections have been applied to the observed in- 
frared fluxes in any of the above analysis. If, indeed, the tori in 
these sources are better described by high factors of dumpiness 
(e.g., Honig et al. 2006) with discrete co-spatial clouds being 
responsible for both emission and obscuration, then such correc- 
tions will be invalid. As such, our tight correlation in Fig. 1 using 
uncorrected IR fluxes seems to be an aposteriori fact support- 
ing (but not proving) the non-applicability of extinction correc- 
tions. Some previous studies have applied extinction corrections 
(Krabbe et al. 2001), while others have not (Lutz et al. 2004). 
In this section, we comment on the extent that canonical cor- 
rections would have on the inferred luminosities and the Lir:Lx 
correlation. 

For our entire sample, the median \ogNn ~ 23.1. Using the 
standard interstellar extinction law (Rieke & Lebofsky 1985), 
the extinction Ai2^ m =0.03Av. A typical Galactic gas:dust ratio 
(Bohlin et al. 1978) then implies a median increase of the L M ir 
values by a factor of 7 for our sample, if one were to correct the 
mid-IR fluxes for such reddening. On the other hand, for sources 
with higher obscuration, this correction will far over-estimate the 
intrinsic power. For instance, the correction factor for a source 
obscured by a gas column of A^h=4.5x10 23 cm -2 (the median 
value of our obscured sample) is about 1000; a corresponding 
increase in Lmir would end up pushing the Bolometric lumi- 
nosities well above the expected Eddington luminosities. 

A possible solution may be to use smaller dust:gas ratios. 
Av/Wh values lower by a factor of ~ 10 than the Galactic value 
have been inferred for nearby Seyferts (Maiolino et al. 2001); 
this would moderate some of these large corrections, but would 
still over-estimate the corrections for the extreme Compton-thick 
cases. We note that most of the sources with the lowest Lmir/^x 
ratios (the ones highest above the correlation line in Fig. 1) are 



all highly obscured and one of these (NGC 7674) is Compton- 
thick. Recently, Zakamska et al. (2008) have found that deep 
Si absorption is more likely to occur in AGN with the highest 
obscuration only. The median Si optical depth for their likely 
Compton-thick sources is T9.7 « 1 .5, but with a significant spread 
of T9.7 ~ 1 in their small sample (similar spreads have been in- 
ferred by Shi et al. 2006 and Hao et al. 2007). Assuming our 
highly-obscured sources also have similar IR opacities, and us- 
ing a standard Milky Way opacity curve (Draine 2003) between 
9.7 and 12-/im, correction factors of AlogLjyiiR=0. 1-0.4 (corre- 
sponding to A v =7-35) should be applied to obtain the intrinsic 
12.3-^m power of Compton-thick sources. Adopting such cor- 
rections has only a small effect on the fitted correlation parame- 
ters in Table 2, changing them by less than their 1-cr uncertain- 
ties. 

To reiterate, we do not apply any IR extinction correction. 

6. Compton-thick sources 

Reliably identifying and gauging the intrinsic powers of 
Compton-thick (CT) AGN remains a difficult issue irrespec- 
tive of source redshift (e.g. Comastri et al. 2007; Nandra & 
Iwasawa 2007; Alexander et al. 2008). In most sources with 
extreme obscuring columns, indirect probes have to be em- 
ployed and allowances made for extinction and contamination 
on extended scales around the nucleus (where the indirect flux 
emerges from). It may well be the case that both the X-ray and 
IR fluxes require additional corrections for obscuration and dust 
extinction. So there is no apriori reason to expect such sources 
to fit the correlation. 

It is therefore encouraging, if somewhat surprising, that our 
eight CT AGN with \ogNu>2A2 do not deviate significantly 
from the best-fit correlation. Due to the typically large uncertain- 
ties in their luminosities (see details in Appendix of this paper 
and of Paper II), CT AGN do not have a substantial effect on 
the correlation fit. But even ignoring these errors, the absolute 
scatter in log(LMiR/i-x) is still moderate (see Fig. 2). 

The importance of our study lies in the fact that the mid- 
IR:X-ray correlation provides not only an independent, but also 
accurate (with small scatter) way to probe AGN, irrespective of 
any obscuration. Simple, high-resolution mid-IR imaging now 
provides a new proxy for AGN power including CT AGN. We 
can test this by searching for high-resolution mid-IR imag- 
ing observations of other CT Seyferts. Recently, Honig et al. 
(2008b) obtained high spatial-resolution iV-band spectroscopy 
with VISIR of a small sample of local AGN and discovered a 
strong mid-IR emission line Baldwin effect. Included in their 
sample are three CT AGN, of which NGC 1068 is already part 
of our sample. The second object, Circinus, is a special case 
(as discussed by these authors) and also showed significant ex- 
tended emission; therefore, we do not consider it here. This 
leaves ESO 428-G014, which was first suggested as being CT 
from BeppoSAX observations by Maiolino et al. (1998). Honig 
et al. (2008b) quote a 12-yum flux of 0.22 Jy => logL M iR=42.5 1 
at a distance of 22.3 Mpc. Our mid-IR:X-ray correlation from 
Eq. 2 then implies logLx=42.39±0.06 - the error being a 1-cr 
Monte Carlo uncertainty. Comparing with recent determinations 
of the intrinsic X-ray power based on the Fe Ka and the opti- 
cal [Om] /15007 emission lines, Levenson et al. (2006) estimate 
logL x =42.40^2.81 for ESO 428-G014 (after applying a small 
distance correction for consistency). This overlaps with the pre- 
dicted logLx from our correlation to within 1-cr. 

The mid-IR and X-ray observations for ESO 438-G014 are 
over-plotted in Fig. 1 and provide additional validation of our 
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correlation in the CT regime. Of course, this result is based on 
a relatively-small sample of eight CT sources only, and further 
observations of an enlarged sample will be important to make 
these conclusions definitive. 

There is one final noteworthy observation regarding star- 
formation around CT AGN. If nuclear starbursts occur in disks 
aligned with the obscuring AGN tori, they should lie along our 
line-of-sight to the core and contaminate the nuclear flux, ir- 
respective of angular-resolution; in this case a systematically 
higher Lir/Lx is expected in our CT point-like cores, which we 
do not observe. This seems to suggest the absence of any on- 
going star-formation (even out to several kpc) around local CT 
sources at the present epoch, at least to within the power limits 
that we discuss in § 4.3. 

7. Summary 

In this work, we suggest that high-resolution mid-infrared con- 
tinuum photometry is an accurate proxy for the intrinsic X-ray 
luminosity of local Seyferts. The emitted powers at ~12-yum and 
2-10 keV are intimately correlated over three orders of magni- 
tude in luminosity, and evidence suggests that this extends into 
the quasar regime. Interestingly, it has been noted in a preprint 
by Grossan (2004, based on high-resolution Keck imaging) that 
such a mid-IR:X-ray correlation may be valid for low-luminosity 
AGN as well. The suitability of the 10-12 (im band for AGN 
selection has been emphasized by other earlier studies as well 
(e.g. Spinoglio & Malkan 1989; Alonso-Herrero et al. 2001). 
Our tight correlation now extends this to the highest spatial res- 
olution available, for Seyferts of all types including very heavily 
obscured sources, and places the conclusion on a precise, quan- 
titative footing with very small scatter. 

Combining our sample with Papers I and II, we have pre- 
sented the most unbiased mid-IR fluxes to date of 41 AGN (not 
including NGC 1068, which was published by others) - 16 of 
these in the current paper. Selecting a sub-sample of sources 
based on a physically-motivated (but empirically-calibrated) 
threshold of physical scale resolution related to the torus size, 
we find the correlation to be especially tight: Lmiro^x 1 ' 11 * 007 
(Eq. 2) with additional systematic variation over the range of 
1.02-1.21. We also find a well-defined range in mid-IR:X-ray 
spectral indices (aix=-1.10±0.01) for this sample, and show 
that mid-IR bolometric corrections (L Bo \/Lmir, indirectly de- 
rived from X-rays) range over ~ 10-30 for the Sy luminosity 
range that we probe. 

It is important to reiterate that even without any selection 
based on resolved scale, our full dataset still has a well-defined 
distribution in Lmir-vs.-Lx (Fig. 2). In fact, if one has access to 
mid-IR imaging of distant AGN (in effect, low-resolution imag- 
ing) and wishes to use this to determine intrinsic X-ray pow- 
ers, then using the parameters for the correlation fit to all 42 
sources (see All' row in Table 2) is probably more appropri- 
ate. Similarly, if one wishes to predict total mid-IR fluxes for 
small-telescope imaging based on known Lx values, the same 
parameters ought to be used. On the other hand, if the goal is a 
determination of the intrinsic mid-IR and X-ray emission of the 
AGN+torus, then the parameters for the well-resolved sample 
(i.e. Eq. 2) can be used directly. 

There are several other implications of our work, and re- 
sultant questions, that we have touched upon only briefly in 
Papers II and herein. That mid-IR nuclear starbursts cannot 
be energetically-dominant within the central ~ 70 pc of local 
Seyfert nuclei was discussed in § 4.3. Our tight dispersion im- 
plies that SBs account for no more than about 40% of the mid-IR 



emission that remains unresolved with the VLT, on average; this 
is to be considered as a strong upper-limit, since other effects 
that increase the dispersion must also be present. The Herschel 2 
satellite (Poglitsch et al. 2006, scheduled for launch in 2009), 
with its large 3.5-m diameter mirror and far-IR detectors, will 
provide an excellent opportunity to test this result with sensitive 
(spatial and spectral) deconvolution of SBs from AGN, better 
than any other space observatory to date. On a more immedi- 
ate timescale, the Infrared Camera (IRC) on Japan's Akari 3 IR 
satellite has several filters sensitive around the 12-yum regime, 
one of which (S9W) has been used for an all-sky survey more 
sensitive than IRAS by an order of magnitude (Murakami et al. 
2007; Onaka et al. 2007; Ishihara et al. 2008). Our even-higher- 
resolution 8-m photometry and Seyfert correlation will give a 
very accurate scaling for removal of contamination in broad- 
band spectral modeling of AGN with these and other future mis- 
sions. 

In this regard, we have not speculated on what the resolved 
scale threshold of -400-1700 

r sub (§ 5.1) could correspond to 
physically. If this is indeed ~10 times larger than typical outer 
torus radii found in recent models (Honig et al. 2006), then it 
may instead correspond to a transition between inner and outer 
nuclear star-forming components (cf. Davies et al. 2007, who 
found a inner star-formation scale of ~ 50 pc). But more work is 
needed to verify this threshold and understand its origin. 

Next, the increase of bolometric corrections - in other words, 
a decrease of Lmir/^BoI - with increasing LboI is qualitatively 
similar to that inferred in the quasar regime by other authors 
(e.g. Maiolino et al. 2007; Treister et al. 2008) and may be inter- 
preted in terms of a decrease of covering factor of dust for more 
powerful sources. Our correlation extends this result at least one 
order of magnitude into the Seyfert regime. Comparison with the 
steeper slope of the X-ray bolometric corrections with luminos- 
ity (Fig. 5) has consequences for the covering factor of obscuring 
gas. For instance, the nearly-linear luminosity correlation may 
imply that the gas as the dust are closely coupled together at all 
powers. But we have not explored this in depth. 

Finally, the fact that all types of Seyferts and Compton-thick 
AGN follow the same correlation suggests a similar origin of the 
mid-IR radiation in all these AGN classes. A non-thermal com- 
ponent is probably not dominant in most cases, as discussed in 
§ 4.2. It may be best explained in the context of clumpy AGN 
tori discussed in Paper II (and references therein). This can also 
be emphasized by plotting the ratio of mid-IR:X-ray luminosi- 
ties against the obscuring columns Nn- See Fig. 8 which is an 
update of Fig. 2 in Paper II. Any difference in the distribution 
of log(LMiR/ix) values for obscured sources with logA^H>22, 
as compared to less obscured ones, is significant only at 78 per 
cent according to a KS test. Furthermore, the corresponding dis- 
tributions for the well-resolved sub-samples (obscured vs. un- 
obscured) are completely indistinguishable. If the gas column 
density is a measure of the torus inclination angle to the line- 
of-sight, then smooth-dust torus models would predict a lower 
log(LMiR/ix) at high A^h values, which we do not see. 

These issues must be accounted for when simulating radia- 
tive transfer through AGN tori. We leave detailed study to future 
work including modeling and observations of enlarged and more 
complete samples. 
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Fig. 8. Log luminosity-ratio as a function of obscuring column density determined from X-ray observations (for all sources excluding ESO 323- 
G032). (r, o>) = (0.24, 0.38) and (0.37, 0.32) for the 24 obscured (log/VH>22) and 17 unobscured sources respectively. The corresponding values 
for the well-resolved sample are (0.17, 0.25) and (0.12, 0.21) with 15 and 7 sources respectively, with the two classes being indistinguishable 
according to a KS test. Nn uncertainties are not included in the plot, but this is irrelevant from the perspective of the absence of any trends. See 

§7. 
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Appendix A: Notes on individual sources 

A.1. Compton-thin sources 

A.1.1. SWIFT J01 38.6-4001 and SWIFT J0601. 9-8636 

The Sy type classification for SWIFT J0138.6-4001 comes 
from Bodaghee et al. (2007, and references therein); that for 
SWIFT J0601 .9-8636 from Landi et al. (2007). The two sources 
were only recently identified as highly obscured AGN with very 
small scattering fractions by Ueda et al. (2007). Theirs is the best 
broad-band analysis (combining Suzaku and Swift data), and we 
use their luminosities directly after correction for transmission 
efficiency. Some variability (by factors of 1.6 and 0.5 respec- 
tively) as compared to the results of Tueller et al. (2008) was 
discussed by Ueda et al. (2007); we thus assign factor of 2 lumi- 
nosity uncertainties. 

SWIFT J060 1.9-8636 was the least significant of all our 
VISIR targets, barely detectable at the 5-cr level. A single source, 
offset by A(Right Ascension)=-l"6, A(Declination)=+2"3 
from the nominal position taken from the NASA Extragalactic 
Database, was found in several filters. Other VLT optical ac- 
quisition images obtained by us as part of a separate program 
confirm that this source coincides with the peak of the optical 
emission from this highly-reddened, edge-on galaxy. 

Both targets were interpreted by Ueda et al. (2007) as be- 
ing buried within geometrically-thick tori, and it is interesting 
to note that one of them (SWIFT J0138.6-4001) lies above the 
best-fit correlation line in Fig. 1 (i.e. is comparatively-weak at 
12-/im, for its X-ray luminosity), possibly implying that hot dust 
clouds are smothered within the thick tori and not observed di- 
rectly. Sensitive observations over a broad wavelength range will 
be required to confirm the nature of this source. 

A.1.2. ESO209-G12 

Sazonov et al. (2007) state Fn-60 kev = 1.66 x 10" 11 erg s" 1 
cm -2 => logLn_6o = 43.75 and logA^n<22 based on their 
INTEGRAL and ROSAT analysis. Using xspec (Arnaud 1996) 
and assuming an unabsorbed power-law (PL) model with V = 1 .7 
and the above normalization also agrees well with L20-100 stated 
by Bassani et al. (2006) when extrapolated to higher energies. 
This PL then implies L2-10 = 4.45 X 10 43 erg s -1 . The model 
straddles the 20^0 keV luminosities stated by Bird et al. (2006) 
and Bird et al. (2007), thus giving a simple uncertainty estimate 
between these. 

A.1.3. NGC 3081 

Beckmann et al. (2007) find, from a combined Swift and 
INTEGRAL analysis, T= 1.8 ±0.2 with L20-100 = 8.1 X 10 42 erg 
s -1 . They also find significant variability of a factor of a few in 
the BAT 14-195 keV lightcurve over a 15 month period. 

Fitting an absorbed PL above a few keV (the softer flux 
is fit by weak scattering or a blackbody) to spectra extracted 
from the Tartarus database gives A^h=68(±5)x 10 22 cm -2 with 
£2-10 = 4 - 6 x 10 42 erg s ; the variation being that between 
the ASCA SISs. Extending the de-absorbed PL to 200 keV gives 
£44-195 = 9 - 14 x 10 42 erg s -1 , well matching the more re- 
cent values quoted by Markwardt et al. (2005) and Ajello et al. 
(2008). This suggests that variability may not be very strong. 

On the other hand, BeppoSAX is known to have observed the 
source in a low flux state (Maiolino et al. 1998). We take the 
lower value of the luminosity range found in the above ASCA 
fits, with a large error of logL=0.3 as a mean value. 



This source was also imaged by Krabbe et al. (2001) at 
1"3 resolution. The full N-band flux quoted by them is consis- 
tent with our narrow-band PAH2 flux (assuming a flat SED in 
AFx units), implying the absence of extended emission on ~ 1" 
scales. 

A.1.4. ESO263-G013 

The only good-quality X-ray spectrum published is the Swift 
XRT analysis by Landi et al. (2007), who find ^h=4+1 xlO 23 
cm -2 absorbing a PL with T « 2.32. The observed flux is F2-10 = 
2x 10~ 12 erg s -1 cm -2 and they quote an intrinsic L2-10 =4 — 
7 x 10 42 erg s -1 (converted to our cosmology). But extending this 
steep PL to higher energies underpredicts other INTEGRAL and 
Swift fluxes by factors of ~ 7 - 10 (Bird et al. 2007; Ajello et al. 
2006). Recent XMM-Newton data (by M. Ajello, priv. comm., 
and some of us) analysed in conjunction with the Swift BAT 
data are instead fit with a harder photon-index T « 1.8, resulting 
in reasonable agreement (to within a factor of ~ 2) with fluxes 
over 2-100 keV and an intrinsic L2-10 ~ 2 x 10 43 erg s -1 ; we 
prefer this latter broad-band analysis. 

A.1.5. NGC 4388 

Combined analysis from several missions shows a simple PL 
above a few keV absorbed by a column N\\-2.1 x 10 23 cm -2 , but 
with highly variable fluxes and obscuring column (Beckmann 
et al. 2004; Cappi et al. 2006; Iwasawa et al. 2003; Risaliti 
et al. 2002). Using the mean of the large range found in XMM- 
Newton by Cappi et al. (2006) and recent INTEGRAL results by 
Beckmann et al. (2006), for a distance of 16.7 Mpc, results in a 
mean logL2-io = 42.24, to which we attach a large 1-cr error in 
dexof0.3. 

A. 1.6. ES0 323-G32 

This source is classified as a Sy 1 .9 in the catalogue of Veron- 
Cetty & Veron (2001). The only reliable published X-ray mea- 
surements happen to be the INTEGRAL observations of Bird 
et al. (2007), who find F 2 o-40 = 6.81 x 10" 12 and F40-100 = 
1.51 x 10" 11 erg s -1 cm -2 . In the absence of other information, 
we parametrize this with a simple PL, T as 1.2, which extrapo- 
lates to F2-10 = 3.7 x 10" 12 erg s" 1 cm" 2 , or L 2 -io = 2.0 x 10 42 
erg s -1 . One could also predict the 2-10 keV power based on 
forbidden optical emission line measurements. Correcting the 
[Om] /15007 emission line flux from Gu et al. (2006) for a 
small extinction factor from the same paper, and using the type- 
1 </ r 2-io/^iomj> relation from Panessa et al. (2006), we get 
L2-10 ~ 3.3 x 10 42 ergs" 1 . 

A.1.7. NGC 4992 

The Sy classification of this source comes from Bodaghee 
et al. (2007) and references therein. Comastri et al. (2007) 
have presented a broad-band Suzaku spectral analysis, and in- 
fer an extreme reflection component with R > 5 and a possi- 
ble relativistically-blurred Fe line, indicative of disc-dominated 
reflection. This can be naturally explained by relativistic light- 
bending on to the inner disc, which preferentially enhances 
reflection and suppresses the incident power-law at infinity 
(Miniutti et al. 2004). The fraction of AGN which are described 
by this light-bending model is still unknown, though increas- 
ing numbers of sources are being discovered and the hard X-ray 
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background spectrum can also accommodate a non-negligible 
fraction of light-bent source in the Universe (see for instance 
Miniutti et al. 2004, Gandhi et al. 2007, Miniutti et al. 2007 and 
references therein). 

Although the majority of intrinsic power-law photons may be 
bent towards the event horizon and never escape the black hole's 
gravitational sphere of influence, the relevant luminosity in our 
case is the X-ray flux that is incident on the distant torus clouds, 
i.e. the (absorption-corrected) power observed at infinity. Using 
the model of Comastri et al. (2007) yields L 2 -io = 1.7 x 10 43 
erg s _1 . If part of the high reflection is due instead to decreased 
transmission efficiency by heavy obscuring material or due to 
strong time variability, then this estimate will be a lower limit 
(cf. Ueda et al. 2007). 

A.1.8. NGC 6300 

This is one of a sample of 'changing-look' AGN whose spec- 
tra exhibit transitions between Compton-thick and Compton-thin 
obscuration on timescales of years (Matt et al. 2003). It was 
found to be reflection-dominated in 1997, but has since been in a 
Compton-thin state (RXTE and BeppoSAX analysis; Guainazzi 
2002). Assuming a Compton-thin column of Nn=2 x 10 23 cm -2 
(from the XMM-Newton analysis of Matsumoto et al. 2004) ob- 
scuring a PL with r = 2.42 (Beckmann et al. 2006) and normal- 
izing to JEM-X and ISGRI 2-100 keV fluxes of the latter pub- 
lication, the implied intrinsic logL2-io = 42.76. But, in addition 
to its changing-look nature, the source also exhibits large flux 
variability (Awaki et al. 2006) and a much lower luminosity of 
logL2-io = 41.8 was derived by Matsumoto et al. (2004) from 
their XMM-Newton analysis. We use a mean logL2-io = 42.3. 

A.1.9. ESO 103-G35 

Molina et al. (2006) analysed the broad-band INTEGRAL spec- 
trum of this source and located an exponential cut-off to the 
X-ray PL at E cut ~ 70 keV, with no substantial flux change 
from previous BeppoSAX observations (Risaliti 2002; Wilkes 
et al. 2001). This source has shown moderate and fast Nh vari- 
ability on time scales of several months (Risaliti et al. 2002). 
Using model realizations in xspec, we find that the PL model 
of Molina et al. (2006) with F = 1.78, JV H =2x 10 23 cm" 1 and 
F2-10 = 3 x 10 -11 cm -2 is consistent with the higher energy 20- 
100 keV flux when extrapolated, but slightly underestimates (by 
~ 30 per cent) the latest 20-40 and 40-100 keV fluxes from Bird 
et al. (2007). On the other hand, the Swift BAT 14-195 keV lu- 
minosity is over-estimated by a factor of a 2. We assume the 
base 2-10 keV luminosity of the Molina et al. model, with an 
1-cr luminosity error of 0.15 in dex. 

A.1.10. NGC 6814 

This source is an unabsorbed Sy 1.5, with logA^n<20.7 
(Reynolds 1997). It is known to be highly variable, even after 
accounting for a neighbouring CV which contaminated early 
observations (e.g. see the results of long-term RXTE monitor- 
ing described by Mukai et al. 2003). Using the PL model 
with T = 2.48 from Beckmann et al. (2006), and extrapolating 
from their INTEGRAL 20-100 fluxes to below 10 keV sug- 
gests that JEM-X should have detected the source, which it did 
not. Instead, some spectral curvature probably needs to be in- 
voked. Including a simple pexrav (Magdziarz & Zdziarski 1995) 
parametrization with r = 1.6, R = 1 and £ cut = 150 keV agrees 



with 20^0 keV and 40-100 keV fluxes of Bird et al. (2007), 
agrees with the mean published RXTE observations mentioned 
above, and also matches the < 10 keV photon-index found from 
ASCA analysis of data in the Tartarus database. This implies a 
mean L2-10 = 1-33 x 10 42 erg s" 1 . A large variation of 0.4 in dex 
accounts for the source variability. 

A.2. Compton-thick sources 
A.2.1. NGC 1068 

NGC 1068 is one of the best-studied of obscured AGN; hence 
our decision to include it on our correlation. In X-rays, it 
is known to be substantially Compton-thick and reflection- 
dominated, with logA^H>25 (Matt et al. 2000; Iwasawa et al. 
1997; Koyama et al. 1989). Correcting from the observed flux 
of X-rays scattered from above the obscuring medium into the 
line-of-sight, or from the flux of the forbidden [Om] A 5007 A 
emission line, intrinsic logL2-io estimates of ~ 42.8 -44 are de- 
rived (Iwasawa et al. 2003; Panessa et al. 2006; Levenson et al. 
2006). 

In the mid-IR, several high-resolution studies have been car- 
ried out from the ground revealing bright, extended structures on 
sub-arcsec scales and beyond (Galliano et al. 2005b; Bock et al. 
2000; Alloin et al. 2000). At the highest spatial sampling, con- 
tinuum core fluxes of 9-10 Jy have been reported around 12-pm 
(Mason et al. 2006; Tomono et al. 2001); these correspond to 
\ogAL A * 43.8. 

A.2.2. NGC 3281 

Vignali & Comastri (2002) identified this source as being 
mildly Compton-thick. Using their model (c) which corrects for 
Compton-scattering in a physically-motivated way => L2-10 = 
1.7 x 10 43 and L20-100 = 2.1 x 10 43 erg s -1 in our cosmology. 
Using this model to infer the 20^1-0 keV and 14-195 keV fluxes 
in order to compare with the recent studies of Bird et al. (2007) 
and Markwardt et al. (2005) suggests logL 2 _io = 43.30 + 0.15. 

Krabbe et al. (2001) detected the presence of extended N- 
band emission beyond their resolution limit of 1" in this source. 
Consistent with this, our Of 35-resolution flux (in units of AF A ) 
is even lower than theirs by « 30 per cent. 

A.2.3. NGC 3393 

Corrections based on the Fe line and on the flux of forbid- 
den optical [Om] emission yield L2-10 ~ 7.2 - 17 x 10 42 erg 
s _1 (Levenson et al. 2006). Combining XMM-Newton analysis 
with previous BeppoSAX data, Guainazzi et al. (2005) find a 
decreased column density (N\\~ 4x 10 24 cm -2 ) with respect to 
the complete Compton-thick suppression (Nn> 10 25 cm -2 ) in- 
ferred by Maiolino et al. (1998). Because only a small part of 
the broad-band X-ray flux is transmitted, correcting the observed 
continuum to estimate the intrinsic luminosity is highly uncer- 
tain given the current data quality on this source. Still, Levenson 
et al. (2006) find reasonable agreement when extrapolating the 
BeppoSAX high-energy data to their Chandra bandpasses with 
r a 2. Suzaku may provide better estimates in the near future. 

A.2.4. IC 3639 

The intrinsic high-energy flux from this source is likely to be 
completely depleted due to heavy Compton-thick obscuration 
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with Nu > 10 25 cm -2 (e.g., Risaliti et al. 1999). Analysis is 
complicated by the fact that most of the < 10 keV X-ray flux ob- 
served is likely related to extended emission (Ghosh et al. 2007; 
Guainazzi et al. 2005). Using instead the [Om] /15007 emission 
line (^[oiii] = 3 .4 x 1 0~ 12 erg s~ 1 cm -2 from Lumsden et al. 200 1 , 
corrected for extinction) as a proxy for the intrinsic power (in 
conjunction with the [Ora]:X-ray relation from Panessa et al. 
2006) gives L2-10 ~ 4x 10 43 erg s _1 . This also matches older 
[Om] measurements quoted by Guainazzi et al. (2005), if we cor- 
rect for extinction using the prescription of Bassani et al. (1999). 

A.2.5. NGC 5728 

From a broad-band Suzaku analysis, Comastri et al. (2007) con- 
clude that this source is Compton-thick with A^h = 2.1 X 10 24 
cm -2 . The unabsorbed flux is F2-10 ~ 5 x 10~ n erg s _1 cm -2 , 
implying L2-10 ~ 8.7 x 10 42 erg s -1 . There is additional low- 
energy spectral complexity due to scattering/reflection, but this 
accounts for only 1-2 per cent of the 2-10 keV power. The 
mildly Compton-thick column suppresses the overall X-ray con- 
tinuum by only ~ 5 - 10 per cent; we include these small correc- 
tions. 

From the raw [Om] /15007 optical emission line measure- 
ment of Gu et al. (2006) along with the type-1 </ r 2-io/^[Oiuj> 
relation from Panessa et al. (2006), we derive L2-10 ~ 1.1 X 10 42 
erg s -1 . Correcting for extinction using the Balmer decrement 
yields L2-10 ~ 3.4 x 10 42 erg s" 1 . 

We use the mean of the above luminosities, which also 
agrees with the extrapolated Swift BAT measurement by 
Markwardt et al. (2005), within errors. 

A.2.6. ESO 138-G1 

Collinge & Brandt (2000) found evidence for Compton-thick 
obscuration in this source, based on their ASCA spectral anal- 
ysis and comparison with [Om] optical emission line flux, but 
the reflection fraction and obscuring column were not well con- 
strained. Suppose we assume that the source PL (F = 2.1) is 
mildly Compton-thick with Nn-2 x 10 24 cm -2 , reflection R=l 
with a typical high energy exponential cut-off at 200 keV. Then, 
normalizing the model to recent INTEGRAL 20-40 keV fluxes 
from Bird et al. (2007) leads to a good match to reported mea- 
surements up to 100 keV (Bassani et al. 2006). In this case, 
logL2-io = 42.9 intrinsically, including mild correction for mul- 
tiple Compton down-scattering (Wilman & Fabian 1999). 

On the other hand, using the Sy 1 relationship between the 
[Om] and 2-10 keV fluxes described by Collinge & Brandt 
(2000), a higher logL2-io = 43.10 is implied. Given the paucity 
of good X-ray data on this source, we assume the mean of the 
above luminosities, with a large uncertainty of 0.3 dex. 



